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Abstract 

In the Wigner framework, one abandons the assumption that the usual canonical commuta- 
tion relations are necessarily valid. Instead, the compatibility of Hamilton's equations and the 
Heisenberg equations are the starting point, and no further assumptions are made about how 
the position and momentum operators commute. Wigner quantization leads to new classes of 
solutions, and representations of Lie superalgebras are needed to describe them. For the n- 
dimensional Wigner harmonic oscillator, solutions are known in terms of the Lie superalgebras 
osp(l|2n) and g[(l|n). For n = 3N, the question arises as to how the angular momentum 
decomposition of representations of these Lie superalgebras is computed. We construct gener- 
ating functions for the angular momentum decomposition of specific series of representations of 
osp(l|6iV) and g[(l|37V), with N = 1 and N = 2. This problem can be completely solved for 
N = 1. However, for N = 2 only some classes of representations allow executable computations. 
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1 Introduction 



The Hamiltonian of the n-dimensional non-isotropic harmonic oscillator with mass m and frequen- 
cies uij (j = 1, . . . , n) is given by 

where the position and momentum operators are denoted by pj and qj respectively. We will look at 
this Hamiltonian in the framework of Wigner quantization, a concept which has been triggered by 
Wigner in p] , but was really introduced much later by Palev [21 S] . Since the Wigner perspective 
has already been considered in detail by Lievens and Van der Jeugt in [5], we will only present their 
results succinctly. For a more thorough deduction of the results, we refer to the aforementioned 
paper. 

Wigner quantization dictates that the canonical commutation relations involving pj and <jj 
should be replaced by less restrictive compatibility conditions (CCs). These arise by imposing the 
equivalence of Hamilton's equations and the Heisenberg equations. For the n-dimensional Wigner 
harmonic oscillator, the CCs are given by 

ih 

[H, qj] = ~ — Pj' i H ,Pj} = ihrnui] q j} 
for j = 1, ... , re. By introducing new operators by 

^• = v^^7m^ (2) 

we can rewrite the Hamiltonian as 

3=1 

In terms of these new operators, the compatibility conditions take the form 

n 

^2ajj[{af,aj},a^]=±2u k af, (3) 
i=i 

for k = 1, . . . , n. Since the position and momentum operators are self-adjoint, we have (a^- )t = aj. 
It turns out that we can find operators subject to the latter hermiticity conditions and to the 
compatibility conditions (|3|) in terms of Lie superalgebra generators. 

There are two known classes of solutions for the non-isotropic case; we can express in terms 
of elements of osp(l|2n) and gl(l|n). For each of these solutions, the spectrum of the Hamiltonian 
in specific Lie superalgebra representations was found in [5]. In their paper the authors give a nice 
overview of the relevant representations and their characters, and they present the energy spectrum 
by means of spectrum generating functions. For a detailed analysis, we refer the reader to that 
paper and the references therein. In this text we will summarize the elements of this paper that 
are useful for our purposes. It should be noted that we will restrict ourselves to the isotropic case 
where u)j = oj for all j = 1, . . . , re. 

When both Lie superalgebra solutions are examined, it is time to move forward to the main 
objective of this paper. The solution of the Wigner quantum system under consideration depends 
on the Lie superalgebra representation V, in which the operators act. So the purpose is to study 
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properties of the Wigner oscillator in different representations, one of which will correspond to the 
canonical case. For the three-dimensional iV-particle Wigner harmonic oscillator, i.e. the case n = 
3N, we want to find the angular momentum/energy contents of Lie superalgebra representations 
of osp(l|6iV) and 0l(l|3iV). In particular, this means we will try to discover which representations 
of so (3) occur in the decomposition of specific Lie superalgebra representations at a given energy 
level. We will use the tool of generating functions to achieve this aim. The aspired results are 
known in the canonical case, and we will compare this case to the new solutions offered by Wigner 
quantization. 

First, we discuss the orthosymplectic case in sections [21 [3] and 01 The osp(l|2n) solution to the 
Wigner problem is discussed in section [2j All of the results in this introductory section stem from [5]. 
In section[3]we explain how the angular momentum contents of the osp(l|6iV) representations V(p) 
can be found by means of generating functions, both theoretically as practically. The actual 
generating functions for representations of osp(l|6) and osp (1 1 12) are computed in section [H We 
also plot the angular momentum/energy contents in so-called (E, j)-diagrams and compare with 
the canonical case in this section. The Lie superalgebra solution g[(l|n) is handled in sections [5l El 
and El In the concluding section, we summarize our main results. 

2 The osp(l|2n) solution 

The orthosymplectic Lie superalgebra osp(l|2n) is generated by its odd elements 6 - (J = 1, . . . , n). 
These paraboson operators are subject to the so-called defining triple relations given by [6] 

[{b% bl}, bf] =(e- OSjibl + (e - (4) 

In these triple relations, j, k and I are elements from the set {1, 2, . . . , n} and ??, £, e E {+, — } (to 
be interpreted as +1 and —1 in the algebraic expressions (e — £) and (e — r/)). The even elements 
of osp(l|2n) are formed by taking anti-commutators {b^,b v k }. 

We can use the paraboson operators to find solutions for the Wigner quantization discussed 
earlier. Indeed, writing as 

aj = bj, a+ = b+, (5) 

with (j = 1,2, ... ,n), we see that the compatibility conditions ((3j) are satisfied using the defining 
triple relations @. The Hamiltonian ([1]) then takes the following form: 

Z 3=1 

In order to obtain that (af)^ = oj, we need to work with suitable representations of osp(l|2n). 
In the paraboson Fock space V^(p) we automatically have (6- )t = bj, which makes this unitary 
irreducible representation of osp(l|2n) an appropriate choice. In [7] the representation V(p) was 
thoroughly investigated, resulting in an explicit basis, matrix elements and character formulas. 
The main theorem of that paper gives the conditions on p for V(p) to be a unitary irreducible 
representation and it states the character of the representation [71 Theorem 7]. 

Theorem 1 The osp(l|2n) representation V(p) with lowest weight (|, . . . , |) is a unirrep if and 
only if p G {1, . . . , n — 1} or p > n — 1. 

The character ofV(p) is given by 

charF (p) = (si ••• x n f' 2 s ^ x ) ( 6 ) 

\,l{\)<\<S>-\ 

The ceiling function [~p] is there to cover the cases where n — 1 < p < n. 
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We have used the notation s\(x) = s\(x±, . . . , x n ) for the symmetric Schur function, which vanishes 
when the length of the partition A, denoted by £(X), exceeds the number of variables n. The length of 
a partition is its number of parts. For a deep introduction to partitions and symmetric polynomials, 
we refer to Macdonald [8]. 

For our purposes, the character formula ([6]) is inadequate. Instead the following equivalent 
formula [7] for p € {1, 2, . . . , n — 1} will be more practical: 

d»V(p) = („■■■ ^ m -_^__ y , (7) 

with 

E (o,P) = XX -1 )^^ 1 '"--'^)- 
In this expression for E(o,p) i the sum is over all partitions of the form 

^ \ai+P 02+P ••• «r+P 

in Frobenius notation, and 

c v = ai + a 2 H h a r + r. 

In the Frobenius notation for partitions, the first and second row denote the lengths of the rows 
and columns in the Young diagram of the partition, counted from the diagonal. For the partition 
rj, a typical shape of the Young diagram is given in Figure [TJ 
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Figure 1: Young diagram of the partition rj in Frobenius notation for r = 3. 



A simple expression for E( 0jP ) exists when p = 1 and p = n — 1. In these cases we have 

charV(l) = ( Xl ■ ■ ■ xn) 1 / 2 — (8) 

and 

char^n - 1) = (x, ■ ..xjW/* — r. (9) 

It is possible to find the spectrum of the Hamiltonian H in the representation V(p). In fact, a 
straightforward technique described in [5] delivers the spectrum generating function spec H , which 
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assigns all eigenvalues of H to a power of t. The multiplicity of an energy level E is then given by 
the coefficient of t . For the osp(l|2n) solution, the spectrum generating function takes the form 

specF = £ £ s x (l,...,l)t^+ k \ (10) 

fc>0 A, |A|=M(A)<I>1 

where |A| is the order of the partition, by which we mean the sum of its parts. The energy levels 
are equidistant and can be written as 

*f = Mf + *), 

with k = 0, 1, 2, The multiplicities of these energy levels are equal to 

M4 P) )= E «a(1, ■■■,!)■ 

A, |A|=M(A)<rPl 

For p = 1 all of the results above coincide with the canonical results. Indeed, p = 1 represents 
canonical quantization. 

3 Angular momentum decomposition of osp(l|2n) 

The main objective of the present paper is to find the angular momentum content of Lie superalge- 
bra representations related to the Wigner quantization of the 3D Wigner harmonic oscillator, both 
for osp(l|2n) and with n = 3N. Both cases are dissimilar with respect to the dimension of 

the representation spaces, so a proper approach is needed to tackle both problems. This asks for a 
small clarification. 

We would like to describe the angular momentum content with a generating function. The 
representation V(p) of osp(l|2ra) is infinite-dimensional, which implies that it is impossible to 
construct a generating function comprising all osp(l[2n) representations. Therefore, our objective in 
the osp(l|2n) case is to construct a generating function for every representation V(p) separately. In 
the g[(l|n) solution, examined from section [5] onwards, the representations will be finite-dimensional. 
In that case, the generating function will contain variables characterizing the gl(l|n) representation. 

We would like to find how the Hilbert space in which the Hamiltonian acts decomposes to so (3) 
representations. In this section we will discuss the osp(l|2ra) case, where n = 3N. In that case we 
can rely on the embedding 

osp(l|6iV) D sp(6iV) D u(3N) D u(3) © u(JV) D so(3) u(l) (11) 

to come up with a generating function that represents the angular momentum decomposition of 
the osp(l|6iV) representation V(p). Some explanation is needed to see why this is the correct 
embedding to use, i.e. to see why the angular momentum operators are elements of this so (3) 
subalgebra of osp(l|6iV). 

3.1 Angular momentum 

For n = 3, our physical system is a three-dimensional harmonic oscillator. In a canonical context, 
the angular momentum operators are defined by M = q x p, or 

3 

M 3 = E e jkiQkPi, (j = 1,2,3), 

k,l=l 
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where ejki is the Levi-Civita symbol. Since the position and momentum operator cannot be assumed 
to commute in Wigner quantization, a logical definition of the angular momentum operators is 



1 3 

Mj = - ^2 e jki{Qk,Pi}, (j = 1,2,3), 
Z U=l 



which, by means of ([2]), can be written as 

M i = ^ E e ^{4,«r}> (j = 1,2,3). (12) 

k,l=l 

The compatibility conditions ([3]) do not contain enough information to lead to commutation rela- 
tions between Mi, Mi and M3. However, in the osp(l|6) solution = b^ one finds 

[M i ,M j ]=ihe ijk M k , (i,j,k = 1,2,3). 

These are the commutation relations of so (3). Now, since the operators Mj are in essence elements 
{b~£,b^} of osp(l|6), they belong to the osp(l|6) subalgebra u(3) from the embedding 

osp(l|6) D sp(6) D u(3) D so(3) © u(l), 

which follows from [3 Proposition 3]. The generalization to osp(l|6./V) is rather straightforward. 
The physical system is now an iV-particle three-dimensional harmonic oscillator. The position and 
momentum operators have a second index a, with a = 1, . . . , N. The angular momentum operators 
of the particle a can be written as 

1 3 -ih 3 

m j> = n E £jki{4k, a ,Pi, a } = -5- E e i^KU' a ;, a }> (i = i> 2 ' 3 )- 
fc,«=l fc,i=l 

The total angular momentum is obtained by adding all the angular momenta of the individual 
particles. Thus we have 

N 

M j = Y, M i^ (j = 1,2,3). (13) 

a=l 

These Mj are elements of the u(3) subalgebra of osp(l|6iV) and satisfy the 50 (3) commutation 
relations. Therefore the angular momentum components generate the so (3) subalgebra of u(3) in 
the following chain of subalgebras: 

osp(l|6iV) D sp(6iV) D u(3N) D u(3) © u(N) D so(3) © u(l). 

The question now is how the osp(l|6iV) representation V(p) decomposes with respect to these 
subalgebras. 



3.2 Decomposing the osp(l|6iV) representation V(p) 

The starting point of the decomposition of V(p) is the character of the Lie superalgebra repre- 
sentation given by equation ([6]). Each Schur- function s\(xi, . . . ,x n ), with A = (Ai, . . . , A n ) is the 
character of an irreducible covariant tensor representation of u(n) [9] and corresponds to the u(n) 
representation with highest weight A, where n = 3iV. In other words, equation ([6|) is a u(3iV) char- 
acter generating function. In contrast, we want the result of our analysis to be a representation 
generating function. In other words, the generating function returns all representations of so (3) 
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that appear in the decomposition of a fixed representation V(p) of osp(l|6iV). By means of an 
example we will try to avoid confusion between both concepts. 

The character of V(p) given by equation ([6]) is, as explained earlier, a u(3iV) character generating 
function. Indeed, it contains the characters of all u(3iV) representations in the decomposition of the 
osp(l|6./V) representation V(p). Such a u(3iV) character however, consist of many superfluous terms 
if one only wishes to know which u(3iV) representations appear. After all, a u(3iV) representation 
is characterized by a partition A, so a term x x = x^ 1 . . . x^ n would suffice instead of s\(xi, . . . , x n ). 

Consider the representation V(2) of osp(l|6) for example. Following equation ([9]) we see that 
the u(3) character generating function takes the form 

xix 2 x 3 (l - X1X2X3) ^ 



(1 - xx)(l - x 2 )(l - x 3 )(l - xix 2 )(l - xix 3 )(l - x 2 x 3 ) ' 

The expansion of this function contains all u(3) characters in the decomposition of the osp(l|6) 
representation V(2). The u(3) character generating function could just as well have been derived 
directly from equation ([6]). The partitions A in this equation have a maximum of two parts, so the 
u(3) representation generating function is created by replacing every Schur function in equation ([6]) 
by its leading term x Xl x 2 2 . We obtain 
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x\x 2 x 3 2^ 2^ x i x 2 = T\ \f\ 7- ( 15 ) 

A^OA^ {l-x 1 ){l-x l x 2 ) 

Every monomial x^x^ 2 in the expansion of this easier looking function corresponds to a u(3) 
representation characterized by the partition A. One can verify that the method described in the 
next paragraph, applied to the u(3) character generating function ()14j) . will indeed give (|15p as a 
result. 

Let us now return to the general case and consider the branching to u(3) @ u(N) in (jlip . The 
substitution 

Xi := UjViz, (j = 1,2,3 and 1 = 1,..., N) (16) 
in equation ([7]) yields a character generating function for u(3) (Bu(N). The factor 

z = t^ 

keeps track of the energy, since the power of z after the substitution (fTUj) in © equals |A| and the 
order of the partition A determines the energy level E { *\ as can be seen from equation ([10]) . We 
will keep using this notation throughout the rest of the paper. 

By now it should be clear that the u(3iV) character generating function ([7]) is not a represen- 
tation generating function. Likewise, after the substitution (fTBI) one does not obtain a u(3) © u(iV) 
representation generating function. Therefore, we need to describe a technique for changing a 
character generating function into a representation generating function. 

From character generating function to representation generating function 

A character generating function for a simple Lie algebra q can generally be written as: 

A 

where the sum runs over a fixed set of integrable highest weights of q, and where each character 
X\{v) is the coefficient of a variable N\ of some sort. Suppose the vector 77 = (rji, . . . , r/ m ) has 
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m components, corresponding to the number of nonzero Dynkin labels of A. The Weyl character 
formula allows us to write the characters as 



J2 e(«0 7] w{X+p) 
rf [[ (1-rj ) 

aeA+ 

where W is the Weyl group, p is the Weyl tool and A + is the set of positive roots of the Lie 
algebra. In order to transform F(rj) into a representation generating function, one has to multiply 
F(rj) by IlaeA + (l — V~ a ) an d keep the terms in the dominant Weyl sector. All of this applied 
to our situation, where the Lie algebra u(n) has the symmetric group as its Weyl group, means 
that we need to maintain the terms in ry A , where A is a partition. One method of doing this, is by 
making the substitution 

C2 c m _i 1 

?7l=Cl77i, Tl2 = —m, ••• r)m-l = Vm-l, Vn 



Cl C m _2 C m _l 



in F(rj) n«eA+(l ~ V °) an d keep all positive powers of c\ y . . . , c m . This comes down to finding 
the term in . . . c m-1 in the power series expansion of the function 



^n (1 - fl (ww) ' 

Several computational software packages have specific methods of finding constant terms in an 
expression. 

For our u(3) © u(N) character generating function, we would have to perform the substitution 
described above for the variables Uj, belonging to u(3), and «/, belonging to u(A r ). Since the next 
step in the decomposition (fTl]) is from u(3) ©u(A^) to so(3) fflu(l), we want the u(N) representation 
labels to be replaced by the dimension of the corresponding u(N) representation in the obtained 
representation generating function for u(3) ©u(A r ). 

Introducing the dimensions of the u(N) representations 

Replacing a term v u = v" 1 . . . v 1 ^ by the dimension of the u(N) representation labeled by v = 
(y\, . . . , un) demands knowledge of a so-called dimension generating function, in which the coef- 
ficient of v v is the dimension of the corresponding u(iV) representation. The dimension of such a 
representation is known [8] and equals s„(l, . . . , 1). Thus, the dimension generating function we 
need is of the form 

5>„(i,...,iK. 

V 

An expression for this u(iV) dimension generating function for general N is not known. However, 
for our purposes the u(2) dimension generating function will be enough. In this N = 2 case we 
have 

Su(xi,x 2 )v u = — — -. (17) 

(1 - XlVi)[l - X 2 Vl)(l - XlX2VlV 2 ) 

Indeed, the summation on the left-hand side can be written as 



00 Ai— A2 



E E E -I 

Ai=A 2 A 2 =0 k=0 



Ai— k Ao+fc Ai A2 
X 2 2 V V 2 ' 
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which simplifies to the right-hand side of (|17p . Therefore, the dimension generating function of 
u(2) is 

(18) 



{l-v 1 f{l-v l v 2 ) 

We now want this u(2) dimension generating function and the previously obtained u(3) © u(2) 
generating function - let us denote this by H2(u,v) - to be "substituted" in one another. Saying 
that two generating functions Fi(X,Y) and F2(X,Z), with common variables X, are substituted 
in each other, means that X x is replaced in either of these generating functions by the coefficient 
of X x in the other. This is achieved by finding the term in X° in either F\(X, Y)F 2 (X^~ 1 ' , Z) or 
F 1 (X(- 1 \Y)F 2 (X,Z), whichever is more easily calculated. Indeed, if F\{X,Y) contains a term 
pi(Y)X x and F 2 (X^ 1 \ Z) includes a term p 2 (Z)X^~ x \ then we find a term p\{Y)p 2 {Z) in the 
product of both functions. 

Here, we have H 2 (u, v) on the one hand and the u(2) dimension generating function Q18p on the 
other hand. Substituting these generating functions in each other is done by taking the constant 
term in the variables v\ and v 2 in the expansion of 

H 2 (u,v) 



(l-"i) 2 (l-ffl 

This replaces the variables v\ and v 2 in H 2 (u,v) by the dimensions of the corresponding u(2) 
representations. Again, finding the constant term can be done by various mathematical software 
packages. 

Angular momentum content 

What is left now is a generating function which models the decomposition of a u(3N) representation 
into u(3) representations. These representations of u(3) are labeled by the variables ui,u 2 ,u^, the 
powers of which represent the partition (Ai,A2,A3) that characterizes the representation. Further 
decomposition to 50 (3) is brought about by the known generating function for su(3) Z> so (3) [10J, 
given here in Dynkin label notation: 

1 + PQJ (19) 



(1-PJ)(1-«J)(1-J«)(1 -<?>)• 



in which J is the so(3) label. Hence, the generating function for u(3) D so(3) in partition notation 
can be written as 

G{ux,u 2 ,u 3 = — — — r~2V 20 

(1 - UlU 2 U 3 ){l - Ui J)(l - U\U 2 J)(l - u()(l - u{u$) 

All factors from equation f)19|) appear in (|20p in accordance with the relation [p, q] = [Ai — A2, A2 — A3], 
except for (1 — u\u 2 u 3 ) in the denominator of (120p . This is explained by the fact that the Dynkin- 
label [p, q] is not influenced when a random integer is added to every part of the partition A. 

Substituting one of these generating functions into the other is done by a similar technique as 
before. One just has to multiply the first generating function, embodying the embedding of u(3) 
in osp(l\6N), by G{u^ , tig , ^3" ) and take the term in u\u 2 v^ in this expression. The resulting 
generating function describes the angular momentum content of the osp(l|6iV) representation V(p). 

4 Generating functions for osp(l|6) and osp(l|12) 

Remember that the goal in the orthosymplectic case is to derive a generating function for each 
representation separately. Such a generating function will be a function of two variables, J and z. 
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The former labels the so (3) content of the representation V(p), while the latter accounts for the 
u(l) part. In fact, each power of z stands for an energy level. 

We can now apply the techniques described in the previous section to derive generating functions 
for the angular momentum decomposition of the representations V(p) of osp(l|6) and osp(l|12). 
However, the osp(l|6) case simplifies drastically as each representation V(p) decomposes to u(3) 
right away. Therefore we can use a different logic to find the desired generating functions. 



4.1 Generating functions for osp(l|6) D so(3) © u(l) 



From the character formula in equation ([6]) one can see that the representation V(p) of osp(l|6) 
decomposes as a direct sum of u(3) representations labeled by a partition A, where A has at most 
three parts. The branching of these u(3) representations can immediately be obtained with the 
help of equation (|20l) . We separate three cases: p = 1, p = 2 and p > 2. 

p = 1 : All partitions in the character formula ([6]) have length 1, so A2 = A3 = 0. It is then 
obvious that the generating function for p = 1 is simply G(z, 0,0), where G(u\, 112,113) is the 
generating function ([6]). Not forgetting the factor z 3p / 2 for the energy we obtain 



7 3/2 



(1 - Zj)(l - Z 2 )' 



(21) 



One can use this generating function to derive the so (3) representations that emerge at energy level 
E^ . This information can be made accessible by means of a table in which the element in row k + 1 
and column j + 1 (counted from the bottom) marks the number of representations J 3 at energy 
level Efp in the angular momentum decomposition of osp(l|6). We call this the (E, j)-diagram of 
osp(l|6) for p = 1. 
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Indeed, the first few terms in the expansion of (|21|) are 



,3/2 



+ J *5/2 + (1 + J2) z m + ( J + J3) z 9/2 + (1 



J 2 + J 4 



z 11 ' 2 + 



We see for example that at energy level Ejp = 9/2Huj, there are two so (3) representations in the 
decomposition of the osp(l|6) representation V(l), characterized by j = 1 and j = 3. Of course, 
these results were already known because p = 1 represents the canonical case. This (E, j)-diagram 
for instance also appears in [11]. 

p = 2 : The partition A now has at most two parts, so A3 = 0. The so(3) u(l) decomposition 
of u(3) representations labeled by such partitions is given by G(z, z, 0). Therefore, we can write 
the generating function for the angular momentum decomposition of osp(l|6) for p = 2 as 



(1 + z 3 J) z 3 



(1 -zJ){l -z 2 J){l -z 2 ){\ -z 4 ) 



(22) 
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As in the previous case, we can generate the (E, j)-diagram of osp(l[6) for p = 2. 
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Let us look at the case Ek = 7Huj, i.e. k = 4 as an example. There are three partitions with two 
parts of order 4, namely (4,0,0), (3,1,0), and (2,2,0). The so(3) representations that emerge in 
these cases can be found by equation (|20p . In total we have 

(1 + J 2 + J 4 ) + (J + J 2 + J 3 ) + (1 + J 2 ). 

This is in accordance with the coefficient of z 7 in equation (|22p . as can be seen from the (E,j)- 
diagram as well. 

p > 2 : Since the length of the partitions in ([6]) cannot exceed the number of variables, we are 
looking in this case at partitions of length at most 3. So the generating function for p > 2 is 

z 3 ^ 2 G(z,z,z), or 

(1 + z 3 J) z 3 P/ 2 



(1 - z.J)(l - z 2 J)(l - z 2 )(l - z 3 )(l - z 4 
The (E, j)-diagram for p > 2 is given by. 



(23) 
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Notice that for the lower energy levels, the cases p = 2 and p > 2 do not differ very much from the 
canonical case. The larger discrepancies are found in higher energy regions. 



4.2 Generating functions for osp(l|12) D so(3) © u(l) 

The previous case might have been very elementary, for osp(l|12) the computations are much 
harder. In fact, it is not practically possible to find generating functions for all osp (1 1 12) represen- 
tation y(p). For V(l) and V(2) however, we are able to follow all the steps from section T3.2I to 
construct the generating function for the angular momentum decomposition. We were unable to 
compute these generating functions for the representations V(p) with p > 3. 



The representation V(l) 

We start with the character of V(l), given by equation ([8]), 
creating a u(3) ® u(2) character generating function: 

(u 2 u 2 u 2 v\V'ifl 2 z ? ' 



and perform the substitution (|16p , thus 



(1 - U\V\Z)(1 - U 2 VIZ)(1 - U 3 Viz)(l - UiV 2 z)(l - U 2 V2Z)(1 - U3V2Z) 
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We need to change this into a representation generating function. To this end, we multiply the 
previous function by 

n n 

where A + and A+ are the positive roots of u(3) and u(2) respectively. Thus we have 

A+ = {(1,-1,0), (1,0,-1), (0,1,-1)} and A' + = {(1, -1)}. 
Therefore we need to multiply our u(3) © u(2) character generating function by 

(1 --)(!--)(! "-)(!--) 

U\ U\ U 2 V\ 



and perform the substitutions 

u\ = au\, u 2 



bu 2 u 3 v 2 
, tt 3 = — , Vi = CVl, v 2 = —. 

a o c 



We want to keep all positive powers of a, b and c, so we multiply our function by 

1 



(1 -a" 1 )(l -6-!)(l -c- 1 ) 

and find the constant term in a, b and c. This is the hardest step to compute. The term in a°6°c° 
factorizes nicely as 

(u\u2u\ VlV 2 ) 3 / 2 Z 3 
(1 — UiV\Z)(l — U\U2V\V 2 Z 1 ) 

This is the u(3) © u(2) representation generating function, in which we want to change the u(2) 
labels v\ and v 2 by the dimensions of the corresponding u(2) representations. Following equation 
(|18p . the u(2) dimension generating function is 



(1 - vi) 2 (l - VlV 2 )' 

Substituting this into the previously obtained u(3) © u(2) representation generating function gives 

UlU 2 Uz z 3 



(1 — tiiz) 2 (l — U\U 2 Z 2 ) 



The angular momentum content is then found by substituting this function and G(ui,u 2 ,u^) from 
equation (I20j) into each other. The resulting angular momentum generating function for the repre- 
sentation V(l) of osp(l|12) is 

(1 + Jz 2 ) z 3 



(l-z 2 ) 3 (l- Jz f 
The (E, j)-diagram shows the result for the first few energy levels. 
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The representation V(2) 

The character of V(2), obtained from ([7]), does not factorize nicely. The angular momentum 
generating function for V(2) can be constructed in precisely the same manner as for V(l). The 
u(3) © u(2) representation generating function in this case equals N/D, with 

,321 33^6 o„, 321„,42. v 6 



N 



1 + u 21 v 21 z 3 + u 221 v 32 z 5 + u 321 v 42 z 6 + u^v M z° - 2u a "v-z" 



u 42 W 
u 43 VV 



u 42 W 



u 422 v 53 z s 



u 431 v 53 z 8 



u 432 v 63 z 9 



2u 532 v f *z w + u 53 W° + u 53 W° + u 63 VV 1 



+ U 64 W 3 + U 85 V°'V 6 



(25a) 



(25b) 



and 

D = (1 - uVz)(l - uW)(l - ^W)(l - i^WXl - u 11 W) 
(1 - U 2 W)(1 - u 21l v 31 z 4 ){\ - u 21 W)(l - u 222 v 33 z 6 ), 

where we have used the notation and similarly for v. The powers of z are 

integers, not partitions. The angular momentum generating function, also quite cumbersome, has 
a numerator equal to 

6 ( l - 2z + 3z 2 + Jz 2 - 2z 3 + AJz 3 + 6z 4 - 7Jz 4 - 3J 2 z 4 - J 3 z 4 - 6z 5 + 4Jz 5 - 2 J 2 z 5 
AJ 3 z 5 + 6z 6 - 8Jz 6 - 3J 2 z 6 + WJ 3 z 6 + 3J 4 z 6 - 2z 7 + WJz 7 + 4JV - 8J 3 z 7 



+ 3z b 



UJz* 



16J 2 / + 13J 3 z 8 



3J 4 z 8 



3J°z 



5 ~ 8 2z 9 + 2Jz 9 + 6J 2 z 9 



+ 6 J 4 z y + 4JV J + z 
+ AJz 11 + 6J 2 z u 



AO 



5Jz 10 + 2JV° + 28JV° + 2jV° - 5JV° + J D z 



6 AO 



12JV 1 +6JV 1 + 2J 5 z 11 



2JV 1 



3Jz 



12 



3JV 2 + 13J 3 z 12 



- 16JV 2 - 14J 5 z 12 + 3J D z 



6 12 



+ iojV 4 



3jV 4 -8jV 4 + 6jV 4 



8 J 3 z 13 + 4J 4 z 13 + 10 J 5 z 13 
- 4JV 5 - 2J 4 z 15 
2 J 6 z 17 + J 5 z 18 + 3J 6 z 18 



2 j6 z 13 +3J 2 z 14 



2 JV 5 + 4jV 5 



- 3J 4 z 16 - 7J 5 z 16 + 6J 6 z 16 + 4J 5 z 17 
while the denominator is 

(1 _ /)4 (1 _ z 2 )2(1 _ z) 2 (1 _ j2z 2 )(1 _ Jz2) 3 (1 _ Jz) 2_ 

Schematically thrown into a (E, j)-diagram this gives 
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for the lower energies. 



5 The £j((l|n) solution 

The elements of the general linear Lie superalgebra gl(l|ra) are denoted by ej k , j, k = 0, 1, . . . ,n, 
where eoj and ejO are the odd generators. The commutation and anti-commutation relations valid 
for this algebra are 

[ey, e kl J = 6 jk e a - (-l) de ^ de ^ 5 u e kr (26) 
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In terms of g[(l|n) generators, solutions of the compatibility conditions ([3]) can be written as 



The hermiticity condition (a^)^ = implies the star condition 

K') f = e j0 . (27) 

The Hamiltonian can then be rewritten as 

i=i 

The unitary representations of gt(l|n) compatible with the star condition (127D are known [12]: aside 
from the typical representations, we have the covariant and contravariant tensor representations. 
Here we are going to work with the covariant tensor representations V\, labeled by a partition A 
with A2 < n. The character of this representation was given by Berele and Regev in [p3]. It is a 
supersymmetric Schur function s\{x\\yi, . . . ,y n ) that can be written as 

charVx = s^x^yi, ...,y n ) 

= ^2 c ^lu s ^(xi)s u >(y 1 ,...,y n ), (28) 

where the coefficients c^ v are the Littlewood- Richardson coefficients and v' is the conjugate partition 
of v, i.e. the partition that is obtained when the Young diagram of v is transposed. The Littlewood- 
Richardson coefficients are integers and arise as coefficients in the expansion of a product of two 
Schur functions as a linear combination of Schur functions. 

Equation (128j) is worth a closer look. First, we note that sJ^c\) vanishes unless the length of 
the partition [i is equal to one, the number of variables of the Schur function. Thus, only partitions 
of the form [i = (r) are allowed. In this case we have 

S( r )(xi) = x\. 

The Littlewood-Richardson coefficients simplify a lot in this case as well. In [8j §5] we find 



c {r)v 



1 if A — v is a horizontal r-strip 
otherwise. 



To explain what a horizontal r-strip is, we must first introduce the notion of a skew diagram. 
Consider two partitions A and v such that Vj < \j for all j. In other words, the Young diagram of 
v is embedded in the Young diagram of A. The set-theoretic difference 9 = A — v is called a skew 
diagram and contains the squares that belong to the Young diagram of A but not to the Young 
diagram of v. If 6 contains at most one block per column, i.e. Q\ < 1, the skew diagram is called 
a horizontal strip. A horizontal strip with r blocks is then called a horizontal r-strip. In Figure [2] 
we find an example of a horizontal 4-strip. 

Combining these results, we find that equation ([28]) can conveniently be rewritten as 



char Vx = s x (xi | y) = ^ x\ ^ s v > (y) , 

r>0 v 

where the second summation runs over all partitions v such that A — v is a horizontal r-strip. The 
first summation is not infinite. Since the Young diagram of the partition A has Ai columns, the 
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Figure 2: The horizontal 4-strip \ — v with A = (5, 4, 2, 2, 1) and v = (5, 2, 2, 1). 



horizontal strip A — v can only have a maximum of Ai parts. So r will never exceed the value Ai. 
Another remark is that r cannot be too small either. If A — v would be an r-strip with r < Ai — n, 
then v would have to be a partition with v\ > n. This would imply that the length of the conjugate 
partition v' is larger than n, which means that the Schur function s u /(y) vanishes. Since \\ — n 
could be negative, we say that r takes values between r* and Ai, where r* is given by 

r * = {° ifAi<n 
[Ai — n otherwise. 

Since A2 < n, horizontal r-strips can always be formed for each r between r* and Ai . In conclusion, 
the character can be simplified further as follows: 

Ai 

charVx = s x {xi\y) = x\^s u i{y), (29) 

where the second sum is again taken over all partitions v such that A — v is a horizontal r-strip. 

Like in the osp(l|2n) case, a spectrum generating function can be produced. We confine our- 
selves to giving the energy levels of the system. They can be written as 

E ( k X) =hu(-^j+r* + k), (fc = 0,...,min(Ai,n)). 

The total number of energy levels depends on Ai and is equal to min(Ai,n) + 1. 

6 Angular momentum decomposition of gl(l|n) 

We will use many of the principles of the orthosymplectic case to find generating functions for 
the angular momentum decomposition of fll(l|n) for n = 3N. One must always bear in mind, 
however, that the g[(l|n) representations V\ are finite-dimensional. Therefore our goal will be to 

A' A' 

create a generating function in which the coefficient of A = A 1 1 . . . A n n represents the angular 
momentum decomposition of V\. We will thus construct a generating function comprising the 
angular momentum decomposition of every 0l(l|n) representation V\. This is different than our 
approach for osp(l|2n), where the generating functions applied to just one representation V(p). 

As before, the angular momentum operators will be part of the so(3) subalgebra of g 1(1 1 3 AT) in 
the chain of subalgebras 

gl(l|3A0 d u(3A0 D u(3) © u{N) D so(3) © u(l) (30) 

However, a little caution is required since it turns out that the angular momentum operators do 
not immediately satisfy the commutation relations of so(3). 
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6.1 Angular momentum 

The angular momentum operators must obviously be denned independently from the choice of the 
Lie superalgebra representation. Thus, just as in equation (fT2j) . we have 

-in 3 

Mj = — £ tjkl{ a k, a Th U = h 2, 3) 
z k,l=i 

for n = 3, i.e. the three-dimensional Wigner harmonic oscillator. In the gl(l|3) solution = eoj 
and aj = ejo one finds 

[M it Mj] = y e ijk M k , k = 1, 2, 3). 

So the operators Lj = 2Mj generate 50 (3). The purpose of angular momentum decomposition is 
mainly finding the spectrum of operators like e.g. M 2 or M3. Obviously, this spectrum only differs 
by a factor from the spectrum of 1? and L3, so finding the 50 (3) content is again a useful problem 
to tackle in this case. 

For n = 3N, the iV-particle three-dimensional Wigner harmonic oscillator, the angular momen- 
tum operators Mj are defined as in equation (jT3j) . Again, apart from a factor 2 these operators 
generate the so(3) subalgebra of g((l|3iV) in the g[(l|3iV) solution a+ = ^/2/(3iV — 1) eoj and 

a J = y?2/(3N — 1) ejQ. Therefore we wish to know how the gl(l|3iV) representation V\ decom- 
poses with respect to the chain of subalgebras (j30|) . 



6.2 Decomposing the g[(l|3iV) representation V\ 

As in the osp(l|6iV) case, the starting point of this decomposition is the character of V\, given by 
equation (|29j) . The Schur-functions s u /(yi, . . . ,V3n) in this character are characters of the u(3iV) 
representations that occur in this decomposition. The branching to u(3) © u(N) representations is 
done by the substitution 

yi:=UjViz, (j = 1,2,3 and I = 1, ... ,N). (31) 

The u(3) and u(N) representations that occur in this branching can be deduced from the following 
known relation: 

S v i (u\Vl, . . . ,U 3 V N ) = ^2g u ', a ,T S a (ui,U2,U 3 )s T (vi, . . . ,V N ), 
cr,T 

which features the Kronecker coefficients g u ',a,r- These are the coefficients that appear when the 
product of two characters of the symmetric group S n are expanded in terms of S n characters: 

a t \ " u' 

^Cp^Cp / j 9u' ,a,T Xq j 

v' 

where a, r, v' and p are partitions of n. Recently King and Welsh [TJj developed a so-called "grand 
generating function" for the Kronecker coefficients. Applied to our context, we can say that when 
u' , a and r are partitions of n, the Kronecker coefficients are the coefficients of the term in z n y u u a v T 
in the expansion of 

All products in this grand generating function run from 1 to the length of the corresponding par- 
titions, which in our case would be 3N, 3 and iV for 2/, a and r respectively. This is already 
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a generating function for the decomposition of u(3iV) to u(3) ®u(iV). But the grand generating 
function is not a u(3) © u(iV) representation generating function since its expansion contains many 
terms of the form y^u^v^ 3 z n in which f^i, fJ>2 and fi^ are not partitions. However, the same tech- 
nique as for the orthosymplectic case will turn the grand generating function into a representation 
generating function. This step is computationally very demanding and can only be performed for 
specific types of u(3A r ) representations. 

Once this step is done, the rest is easy. The u(N) labels v%,...,vn need to be replaced by 
the dimensions of their corresponding u(N) representations as before. The resulting generating 
function then needs to substituted into the generating function ([20]) for the branching u(3) Z> so (3). 

Generating functions for the angular momentum decomposition of the $j[(l|3iV) representations 
V\, where Ai = 1, have been constructed for general values of N in [15] using a different group 
theoretical method. We were able to extend these results to other forms of A, but only for gl(l|3) 
and g[(l|6). For other values of N the computations prove to be too hard. 

7 Generating functions for g[(l|3) and gl(l|6) 

Since the length of the partition A is arbitrary, a generating function in which the coefficient of 
A x = A^ 1 ... is the angular momentum decomposition of the gl(l|3iV) representation V\ would 
have an infinite amount of variables A{. Therefore we choose this angular momentum decomposition 

A' A' 

to be accompanied by A-^ 1 . . . A n n , thus creating a generating function with n = 3N variables. This 
is possible because only the values of A' l5 . . . , A^ affect the angular momentum decomposition of V\, 
as can be seen from the character formula (129p . 

In the previous section we explained how a generating function for the decomposition of a u(3N) 
representation in accordance with the chain of subalgebras 

u(3JV) D u(3) © u(JV) D so(3) © u(l) 

can be created. Let us denote this generating function by H(J, A±, . . . , A n ). The angular mo- 
mentum decomposition of the g[(l|3iV) representation V\ will then be described by the following 
generating function: 

min(Ai,n) 

H(J,A u ...,A n ) z^ +r * J] + M (33) 

i=l 

To see why this is true, we first note that the bottom energy level equals + r*, which explains 
the power of z in (|33p . Each value of r = r* + k then defines a new energy level and is thus 
responsible for an extra factor z. For reasons of clarity, the rest of the analysis will be done for a 
typical partition A, with Ai > n. The same ideas can be adopted in the atypical cases. 

For r = Ai — n, i.e. on the ground energy level, there is only one partition v for which A — v 
is a horizontal r-strip. Its conjugate can be written as v' = (A^, . . . , A^). The corresponding repre- 
sentation of u(3iV) decomposes to so(3) as described by the generating function H(J, A\, . . . , A n ), 
so there must be at least one term 

H(J,A 1 ,...,A n ) 

in the generating function we are trying to describe. In general there are n partitions v such that 
A — v is a horizontal (Ai — n + l)-strip. Their conjugate partitions are of the form 

V = 0\l> • • • ) \ ~ 1) • • • 5 XJi 

where i = 1, . . . , n. The angular momentum decomposition of the u(3iV) representations character- 

A' A' 

ized by these partitions will be in the coefficient of A^ . . . A n n in Ai H(J, A±, . . . , A n ). Therefore, 
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our generating function must also contain the term 

H{J,A 1 ,...,A n ) z{A x + ■■■ + A n ). 

Note that it is possible that A^ = A^ +1 , in which case v' would not be a partition. However, in this 
case the coefficient of A v in H(J, At, . . . , A n ) will be zero. In other words, this "non-partition" 
will not be counted at all. 

The factor {A\ + • • • + A n ) in the term for r = Ai — n + 1 is in fact the elementary symmetric 
polynomial ei(Ai, . . . , A n ). The elementary symmetric polynomial ej in n variables is defined by 
the sum of all possible products of i out of the n variables. We refer to Macdonald [8] for a more 
elaborate introduction to elementary symmetric polynomials. For r = Ai — n + 2 there are typically 
n(n — l)/2 partitions v for which A — v is a horizontal r-strip. An analogous reasoning as before 
shows that they are responsible for a term 

H(J,A 1 ,...,A n ) z 2 A i A i 

l<i<j<n 

in our generating function. This term contains the elementary symmetric function e2(A\, . . . ,A n ). 
Taking all values of r into account, it is not so hard to see that the generating function for the 
angular momentum decomposition of the g[(l|3iV) representation V\ equals 

\x\ n 

H(J, A u ... , A n ) z — +Xl - n J2 ei(A h . . . , A n )z l 

in the typical case where Ai > n. In this expression we find back the generating function for 
elementary symmetric functions, which can be rewritten (see for example Macdonald [8]) as 

n n 

Y,e i (A l ,...,A n )z i = l[(l + A i z). 

i=0 i=l 

For the atypical cases, where Ai < n, we can build up a similar analysis to eventually obtain the 
generating function in equation (f33l) . 

Clearly, the most important part of our problem is finding the u(3iV) D so (3) generating 
function H(J,Ai, . . . ,A n ). However, for g[(l|3) this step is trivial, so the results can be written 
down immediately. 

7.1 Generating functions for gl(l|3) D so(3) 

For the main part, the decomposition of the representation V\ of gl(l|3) following the branching 

gl(l|3) D u(3) Dso(3) 

is described by the generating function G(Ax, A%, A3) given by equation ([20]) . where A%, A2 and 
A3 label the first three parts of the conjugate partition of A. The rest of the generating function 
follows from the previous discussion, equation (i33|) in particular, and depends on Ai. For Ai > 3 
we find 

H +r * (1 + A lZ )(l + A 2 z){l + A 3 z)(l + A\A 2 J) 

Z2 (l-A^A^l-AxJ^l-A^JXl-ADil-AjA 2 )- [ > 

The cases where ^(A') = 2 and £(X') = 1 are easily deduced from this equation by setting A3 = 
and A2 = A3 = respectively. For i(X') = 1 we find back the results from King, Stoilova and Van 
der Jeugt in [15]. 
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The generating function ([34")) allows us to construct (E, j)-diagrams for any fjl(l|3) represen- 
tation V\. Some examples for the typical case are given below, with A a = (3,1,0) and = 
(3,2,2,1,1). 
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The information of these (E, j)-diagrams can be obtained from the coefficient of A v ' a and A v b in the 
expansion of (|34p . for v' a = (2, 1, 1) and v' h = (5, 3, 1) respectively. As a primary difference with the 
canonical case (and the osp(l|6A r ) case in general) we see that there is a finite amount of energy 
levels. Also, at the bottom energy level we see more than one so(3)-multiplet. The first (E,j)- 
diagram, where A = (3,1,0) shows an exception to this remark. We also note that we still have 
equidistant energy levels, and that there are again higher multiplicities of so (3) representations. 

In the atypical cases, the number of energy levels decreases as the length of A becomes smaller. 
For A = (2, 1, 1) we have 
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There are only three energy levels in this case, and we observe that the vertical symmetry of the 
(E, j)-diagram is now gone. 



7.2 Generating functions for gl(l|6) D so (3) 

The case fll(l|3) is deceivingly simple compared to g[(l|6). In fact, it will no longer be possible to 
construct the generating function H(J,Ai, . . . ,A$) for all representations of u(6). More precisely, 
we will only be able to handle the cases £(X') = 1 and ^(A') = 2 completely, where A is the partition 
that characterizes the fll(l|6) representation V\. 



The representation V\ with = 1 

We follow the procedure described in section [6721 to obtain a generating function H(J,A\) which 
will describe the angular momentum generating function for the u(6) representation characterized 
by a partition u', where v' contains of one part only. The starting point is the grand generating 
function (1521 for iV = 2: 



(1 - - s*a-)(l - Ja)(l 



V2 ' 
VI ■ 



(1 - uiviAi)(l - u 2 viAi)(l - u 3 viAi)(l - u\v 2 Ai){l - ti 2 f2^i)(l - u^v 2 Ai) ' 

where we have taken into account the fact that Ai = ^(A') = 1 and therefore only one parameter 
Ai is necessary to describe the u(6) representation. We recognize this function from the V(l) 
representation of osp(l|12), where the starting function was the same. Thus, the rest of the analysis 
can be adopted from that case and eventually we find 



(i-Af)3(i- jA x y 
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The generating function for the angular momentum decomposition of the g[(l|6) representation V\, 
with ^(A') = 1, then follows from equation ([33]) and is equal to 

j*l (l + A lZ )(l + JAl) 

2 (i_ A f)3 (i-jA t y {6b) 

This is again confirmed by the results in |15| . 



The representation V\ with ^(A') = 2 

Compared to the case where i(X') = 1, the denominator of the grand generating function will 
have six extra factors containing A2, and the numerator must have an extra factor (1 — A2/A1). 
From that point on, essentially all computations run along previously traveled paths. Yet, the 
computer has a much harder time performing these computations, and we have the end result (see 
Appendix A) as a witness. Some intermediate results are however interesting. The u(3) © u(2) 
representation generating function in this case equals N/D with 



N = 1 + u 21 v 21 A 21 + u 221 v 32 A 32 + u 32l v i2 A 33 + u 321 v 33 A 42 - 2u 321 v 42 A 42 

- u 421 v 52 A 43 - u 421 v 43 A 52 - u 422 v 53 A 53 - u 431 v 53 A 53 - u 432 v 63 A 54 

- u 432 v 54 A 63 - 2u 532 v m A 64 + u 532 v 73 A 64 + u 532 v M A 73 + u 632 v u A 74 

+ U 643 V 85 A 85 + ^53^10,6^10,6 

and 

D = (1 - u l v l A l ){l - u 2 v n A n )(l - u ll v 2 A ll ){l - u^v^A 2 )^ - u lu v 21 A 21 ) 
(1 - u 22 v 22 A 22 )(l - u 2U v 31 A 22 )(l - u 211 v 22 A 31 ){\ - u 222 v 33 A 33 ), 



(36a) 



(36b) 



where we have used the notation u a = u^u^u^ 3 , and similarly for v and A. This generating 
function was first obtained by Patera and Sharp [16] as a plethysm generating function for two- 
rowed representations of SU(n). In this paper, we have already computed N/D given by equation 
(f36|) independently, not using the grand generating function of King and Welsh. Indeed, this 
generating function and the function N/D, with numerator and denominator given by equation 
(|25p . are similar. The only difference is that A x in (j36|) has been changed into i n equation 
(i25l) . The functions N/D given by equations (f25l) and (j36j) represent the u(3) © u(2) branching of 
all u(6) representations occurring in the representation V(2) of osp(l|12) and V\ (with -^(A') = 2) 
of |6) respectively. The u(6) representations that occur in both cases are characterized by a 
partition with a maximum length of two, so both generating functions must be equal. 

Introducing the dimensions of the u(2) representations in our u(3)©u(2) representation generat- 
ing function, and then substituting the result into G(u\, 112,113) gives us the function H( J, A\, A2). 
The unappealing sight of this function forces us to relocate its full expression to Appendix A. We 
can still write our angular momentum content generating function as 

z 2 ^ (l + A 1Z ){1 + A 2 z) H{J, A 1 ,A 2 ). 



Note that for A2 = we must find back equation (|35p . which is indeed the case. It is interesting 
to see what happens when the powers of A\ and A2 are equal. This means that we are looking 
at u(6) representations with character s,/^!,^), where v' is a partition for which both rows 
have equal length. These representations are of interest in complexity theory and in the study of 
qubits [171 [E] • The generating function in this case can be computed out of the previous one by 
making the substitution 

M 

A x = a A x , A 2 = — 

a 
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and then looking for the constant term in a. We find that the u(6) D u(3) © u(2) branching is 
represented by the generating function 

(1 + u 32l v 42 A 33 ) 

(1 - u 2 v ll A ll ){l - u ll v 2 A ll )(l - u 22 v 22 A 22 ){l - u 211 v 3l A 22 )(l - u 222 v 33 A 33 ) ' 

a result which was obtained recently by King and Welsh [14]. This confirms earlier observations 
concerning inner products of Schur functions pjJJ [20] . The generating function Hy =y ( J, A%, A2) = 
N/D, with 

N = 1 — A + 2A 22 + 3JA 22 + 3J 2 A 22 - A 33 - 3JA 33 - 4J 2 A 33 - 3J 3 A 33 + A 44 

+ 3JA 44 - 2J 2 A 44 + 3 J 3 A 44 + J 4 A 44 - 3JA 55 - 4J 2 A 55 - 3J 3 A 55 - J 4 A 55 

+ 3J 2 A m + 3j3j4 66 + 2J 4 A 66 _ j± A 77 + jA^S 

and 

D = (I — A 11 ) 2 (1 - J.4 11 ) 3 (1 - J 2 A 11 ) (1 - A 22 ) 4 , 

describes the further branching of a two-rowed 11(6) representation to so (3). Note that we cannot 
use this generating function to describe the angular momentum decomposition of the g[(l|6) rep- 
resentation V\, for which A is a partition with two columns of equal length in the Young diagram. 
Indeed, solving such a problem would require the angular momentum decomposition of all u(6) rep- 
resentations characterized by a partition u', such that A — v is a horizontal r-strip, with r = 0, 1, 2. 
For r = 1, the partition u' does not consist of two rows of equal length (in fact, v' x = v' 2 + 1), thus 
the generating function Hy i= \' 2 (J, A\, A2) is of no use for u(6) representations characterized by this 
particular partition v' . 

The representation V\ with = 3 

Computationally, this case can only be worked out when A' x = A2 = A3. A three-rowed u(6) 
representation where all rows are of equal length decomposes to u(3) © u(2) in accordance with the 
generating function of the form N/D, with 

N = 1 — u in v 21 A lu - u 32l v 33 A 222 + u 222 v 42 A 222 + u 321 v 42 A 222 + 2u 432 v 54 A 333 
+ u 531 v 54 A 333 - u 432 v 63 A 333 - v 75 u 543 A 444 + u 642 v 66 A 444 - 2u 642 v 75 A 444 

_ ^753^87^555 _ ^852^87^555 + ^753^96^555 + ^963^10,8^666 _ ^10,7,4^12,9^777 

and 

D = (1 - u lu v 21 A lu )(l - u 21 v 21 A lu )(l - u ul v 3 A lu )(l - u 32l v 33 A 222 ) 
x (1 - u 33 v 33 A 222 ){\ - u 41l v 33 A 222 ){\ - u 444 v m A 444 ) 

The same representation of u(6) decomposes to so (3) in accordance with the generating function 
H K = X , 2=X , 3 (J,A 1 ,A 2 ,A 3 ) = N/D, with 

N = 1 + 2JA 222 + 3J 2 A 222 + J 3 A 222 + 2 J A 333 - 2J 2 A 333 - 2J 3 A 333 - 2J 4 A 333 

+ ^444 _ 4j2j4 444 _ 4j3j4 444 + j5 ^444 _ 2 ^555 _ ^2^555 _ ^3^555 

+ 2J 4 A 555 + J 2 A 6m + 3J 3 ^ 666 + 2J 4 A 6m + J 5 A 8S8 

and 

D = (l- A 111 ) 4 (1 - J A 111 ) 2 (1 - J 2 A 111 ) 2 (1 - A 222 ) 3 . 

The argumentation given in the previous section implies again that this generating function does not 
contain sufficient information for the angular momentum decomposition of the g((l|6) representation 
V x , with Ai = \' 2 = A' 3 . 



21 



8 Conclusions 



For a 3D iV-particle Wigner harmonic oscillator, operator solutions exist in terms of generators of 
the Lie superalgebras osp(l|6A r ) and gl(l|3iV). These operators act in representation spaces of these 
Lie superalgebras. Our goal was to find the angular momentum/energy contents of the representa- 
tions V(p) of osp(l|6iV) and V\ of gt(l|3JV). For N = 1, we have been able to construct generating 
functions representing the angular momentum decomposition for all of these representations. For 
N = 2, the computer allowed us to construct only partial results. For osp(l|12) we have generating 
functions for the representations V(l) and V(2), but for other representations the results proved 
computationally too hard. In the gt(l|6) case, we had to restrict ourselves to representations V\ 
for which A had a maximum of two columns. 

By means of the obtained generating functions, we were able to plot the angular momen- 
tum/energy contents in so-called (E, j)-diagrams. These are tables showing the multiplicities of all 
angular momentum values at each energy level. These tables are a practical tool to compare the 
results of the new Wigner solutions to the well-known canonical case. 

For the 1-dimensional Wigner harmonic oscillator, investigated by Wigner in [I], the energy 
levels in the non-canonical solutions are shifted in height but remain equidistant. The (E,j)- 
diagrams in the osp(l|6A r ) solution show that for all representations V(p) we have a similar behavior 
for the angular momentum contents. Apart from the shifted energy levels, the structure of the 
[E, j')-diagrams in the non-canonical solutions is the same as that for the representation V(l) of 
osp(l|6). The main difference is that the multiplicities of the angular momentum representations 
can be higher than 1, a feature not observed in the canonical case. 

In the 0l(l|3iV) solution, the situation is drastically different due to the finite-dimensional nature 
of the representations. There is a finite amount of energy levels and the number of so(3)-multiplets, 
the angular momentum contents, does not increase when the energy gets higher. On the contrary, 
for higher (and lower) energy levels we see less so(3)-multiplets than in the bulk of the energy 
spectrum. 
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Appendix A. The function H(J, Ai, A 2 ) 

The generating function H(J, A±, A2) for the angular momentum decomposition of a two-rowed 
representation of u(6), has the form N/D, with 

N = 1 - A 11 + J A 2 + 2A 21 + 6 J A 21 + 2A 22 + 3 J A 22 + 3 J 2 A 22 + 3.4 31 - 2JA 31 

- QJ 2 A 31 - J 3 A 31 + 2A 32 - QJA 32 - 8J 2 A 32 - 6J 3 A 32 - A 33 - 3JA 33 

- 4J 2 A 33 - 3J 3 A 33 - 6JA 41 - 3A 42 - 6JA 42 - 6J 2 A 42 + 5 J 3 A 42 + 3J 4 ^ 42 

- 2 J A 43 + 4 J 2 A 43 + 8 J 3 ^ 43 + 6 J 4 A 43 + ^ 44 + 3 J A 44 - 2 J 2 A 44 + 3J 3 A 44 
+ J 4 A 44 - A 51 + 3 J 2 A 51 + 4JA 52 + 6J 3 A 52 + 3A 53 + 2 J^ 53 + J 2 A 53 

+ 3 J 3 A 53 - 4J 4 A 53 - 3J 5 A 53 - 2JA 54 - 8J 2 A 54 - 6J 4 A 54 - 2J 5 A 54 

- 3JA 55 - 4J 2 A 55 - 3J 3 A 55 - J 4 A 55 + 2 J A 61 + A 62 + J 2 A 62 - 2A 63 

- 10JA 63 - 10J 2 A 63 - 3A 64 - 11JA 64 + 4J 2 A 64 + 5J 3 A 64 + 6J 4 A 64 + 3J 5 A 64 
+ J 6 A 64 + 4 J 2 A 65 + 6 J 3 A 65 + 8 J 4 A 65 + 2 J 5 A 65 + 3 J 2 A m + 3 J 3 ^4 66 

+ 2J 4 A m - J 2 A 71 - 2JA 72 - 2J 3 A 72 - 2A 73 + J A 73 + 8J 2 A 73 + 8J 3 A 73 

+ 6 J 4 A 73 - 2A 74 + 8 J A 74 + 30J 2 A 74 + 18J 3 A 74 + 8 J 4 A 74 - 4J 5 A 74 + 8 J A 75 

+ 15, J 2 A 75 + 12J 3 A 75 + 3 J 4 ,4 75 - 7J 5 A 75 - J 6 A 75 - 2J 3 A 76 - 4J 4 A 76 

_ 4J 5 A 76 _ j4 A 77 + j2 ^82 + ^83 + 4J 2^83 + 2J 3^83 + ^84 + ^84 

- 3 J 2 A 84 - 12 J 3 A 84 - 15 J 4 A 84 - 8J 5 A 84 + AJA 85 - 8J 2 A 85 - 18 J 3 A 85 
_ 30J 4 A 85 - 8J 5 A 85 + 2J 6 A 85 - QJ 2 A m - 8J 3 A m - 8J 4 A m - J 5 A 86 

+ 2J 6 A m + 2J 3 A 87 + 2J 5 ^ 87 + J 4 A 88 - 2J 2 A 93 - 3J 3 A 93 - 3J 4 A 93 

- 2JA 94 - 8J 2 A 94 - 6J 3 A 94 - AJ 4 A 94 - A 95 - 3JA 95 - 6J 2 A 95 - 5J 3 A 95 

- AJ 4 A 95 + 11J 5 A 95 + 3J 6 A 95 + mJ 4 A m + 10J 5 ^ 96 + 2,fA m - J 4 A 97 

- J 6 A 97 - 2J 5 A 98 + J 2 A 10,4 + 3 J 3 A 10 ' 4 + AJ 4 A 10 ' 4 + 3 J 5 A 10 ' 4 + 2 J A 10 ' 5 

+ qj 2 A 10,5 + 8 J 4 A 10 ' 5 + 2J 5 A 10,5 + 3 J A 10 ' 6 + 4 J 2 A 10 ' 6 - 3 J 3 A 10 ' 6 - J 4 A 10 ' 6 

- 2J 5 A 10,6 - 3J 6 A 10 ' 6 - 6 J 3 A 10 ' 7 - AJ 5 A 10 ' 7 - 3 J 4 A 10 ' 8 + J 6 A 10,8 - J 2 A 11 ' 5 

- 3J 3 A 11,5 + 2 J 4 A 11 ' 5 - 3J 5 A 11,5 - J 6 A 11,5 - 6J 2 A U ' 6 - 8 J 3 A 11 ' 6 
_ 4 J 4 A n ' 6 + 2 J 5 A 11 ' 6 - 3J 2 A 11,7 - 5 J 3 A 11 ' 7 + 6J 4 A 11,7 + 6 J 5 A 11 ' 7 
+ 3 J 6 A 11 ' 7 + 6 J 5 ^ 11 ' 8 + 3 J 3 A 12 ' 6 + 4 J 4 ^ 12 ' 6 + 3J 5 A 12,6 + J 6 A 12,6 
+ 6J 3 A 12 ' 7 + 8J 4 ^ 12 ' 7 + 6J 5 A 12 ' 7 - 2 J 6 A 12 ' 7 + J 3 A 12 ' 8 + 6J 4 A 12 < 8 
+ 2J 5 A 12 < 8 - 3J 6 ^ 12 ' 8 - 3J 4 A 13 < 7 - 3J 5 ^ 13 ' 7 - 2J 6 A 13 ' 7 - 6J 5 ^ 13 ' 8 

- 2J 6 A 13 ' 8 - J 5 A 13,9 + J 6 ,4 14 ' 8 - J 6 ^ 15,9 

and 

D = (1 - A ll ) 2 {l - A 2 f{l - A 22 ) 4 {1 - JA) 2 {1 - J^ n ) 3 (l - J 2 A 11 ). 
We have used the notation A v> to denote A^-A^ 2 for practical reasons. 
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